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Manuel Valfn-Rodriguez, Antonio Puente, and Llorenc Serra 
Departament de Fisica, Universitat de les Illes Balears, E-07071 Palma de Mallorca, Spain 



(N 
O 
O 
(N 

C 
3 



(N 



c3 



-a 
o 

<N 
> 

oo 
m 

o 

(N 
O 

-i— > 
c3 



i 

c 

o 
o 



X 



Enrico Lipparini 

Dipartimento di Fisica, Universita di Trento, and INFM sezione di Trento, 1-38050 Povo, Italy 

(March 28, 2002) 



The spin-orbit coupling influences the total spin of semi- 
conductor quantum dots. We analyze the theoretical pre- 
diction for the combined effects of spin-orbit coupling, weak 
vertical magnetic fields and deformation of the dot. Our re- 
sults allow the characterization of the quantum dots as spin 
switches, controllable with electric gates. 

PACS 73.21.La, 73.21.-b 



A novel technology based on the use of the electron 
spin, as opposed to the more traditional use of the elec- 
tron charge, is emerging under the name of spintronics. 
Several spin-based electronic devices have already proved 
their importance, even at a commercial level, such as the 
spin- valve read heads. A review of the status of this in- 
cipient field can be found in Ref. Q. 

The spin-orbit (SO) coupling is an essential mechanism 
for most spintronic devices, since it links the spin and the 
charge dynamics, opening the possibility of spin control 
through electric fields Indeed, recent experimental 
and theoretical investigations have shown that the SO 
coupling affects the charge transport and, more specif- 
ically, the conductance fluctuations of chaotic quantum 
dots in a parallel magnetic field 0-0]. It also affects the 
dot far-infrared absorption, introducing peculiar corre- 
lations between the charge and spin oscillating densities 
||. In this work we analyze the combined effects of SO 
coupling, weak vertical magnetic (B) fields, and spatial 
deformation in fixing the spin and other ground state 
properties of model semiconductor dots. 

It will be shown that a sufficiently strong SO cou- 
pling can lead to spin inversion, with an alternating in- 
dependence, similar to the observations from capacitance 
spectroscopy experiments of both vertical J?j and lateral 
quantum dots [0 . Since SO coupling is also active at low 
B's, this mechanism influences the weak field regime and 
can provide an alternative interpretation to the one of 
Ciorga et al. || based on calculations for high magnetic 
fields. 

We consider the SO coupling terms as reviewed by 
Voskoboynikov et al. |J. Assuming a 2D system and 
the effective Hamiltonian formalism, the relevant Dres- 
selhaus contribution for the standard (001) plane of GaAs 
reads 



N 



(1) 



where the a's are the Pauli matrices and P = — ihW + ^A 
represents the canonical momentum given in terms of 
the vector potential A p0|| . The Dresselhaus parame- 
ter Ad is determined by the dot vertical width zq as || 
A D w 7(7r/zo) 2 , with a material-specific constant 7 that 
for GaAs is 7 = 27.5 eVA 3 Q. A SO coupling of Rashba 
type H was also considered, although its contribution to 
the results for GaAs shown below turned out to be neg- 
ligible. 

Neglecting for the moment the Coulomb interaction 
the complete Hamiltonian reads H = Ho + Hd + Hz, 
where Ho consists of the kinetic and confinement ener- 
gies, i.e., 
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Note that the assumed anisotropic confinement will per- 
mit the modelling of dots with varying elliptical shapes. 
The Zeeman term Hz depends on the total vertical spin 
S z , the Bohr magneton [is and the effective gyromag- 
netic factor g*, which for bulk GaAs is —0.44. Namely, 
H z = g*HBBS z . 

Assuming Ho S> Hd 3> Hz and expanding in pow- 
ers of Ad an analytic diagonalization to 0(A^,) in spin 
space is possible with a unitary transformation || H — 
U^HUi, giving the transformed Hamiltonian 
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where we have defined the canonical angular momentum 
operator L z — xP y — yP x - Despite the spin diagonaliza- 
tion, the x and y degrees of freedom in Eq. (||) are still 
coupled through the vector potential in the kinetic en- 
ergy and in L z . With a second transformation each spin 
component can be recast in a separable form using the 
methods of Meyer, Kucar and Cederbaum Specifi- 
cally, defining Hr) = U^jfi^U^n, with rj ==f, J., we obtain 
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where s v = ±1 for r/ =T,4- Assuming, without loss of 
generality, u x > w y the masses and frequencies of the 
decoupled oscillators are 
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with the upper (lower) sign in ± corresponding to k = 
1(2). We have defined in Eq. (|) a spin dependent cy- 
clotron frequency including the SO correction 
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The solution to Eq. (|J) is given by products of x and y 
harmonic oscillator functions which, when transformed 
back to the laboratory frame, yield the desired solutions 
to the original Hamiltonian of Eq. (||). The eigenvalues 
for each spin can be labelled by the number of quanta in 
the x and y oscillators: 
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Figure 1 displays the magnetic field evolution of the 
single-particle energies (0) for a circular dot, except for a 
constant representing the charging energy. These values 
give the chemical potential of the dot for varying elec- 
tron number, as measured for instance in capacitance 
spectroscopy experiments Actually the parabola 

coefficient has been taken from a fit to the experiments 
fl3|| . The spin is indicated in Fig. 1 with the light and 
dark gray tones. In the absence of SO coupling each line 
corresponds to a given spin, except for a very small fluc- 
tuation due to the Zccman energy at some cusps and val- 
leys. In this case the traces arrange themselves in parallel 
pairs of up and down spin. As shown in the two lower 
panels of Fig. 1, the SO coupling produces sizeable up 
and down fluctuations of the spin. These spin inversions 
are due to the level rearrangements embodied in Eq. |t]. 
For z = 100 A, i.e., weak SO coupling, the fluctuations 
start at low magnetic fields and they extend up to B w 
1 T. An even stronger SO (z — 60 A) produces spin 



inversions up to the last level crossing, which marks the 
filling factor v — 2 line Besides, in the latter case 

the traces are no longer paired but, instead, anticorre- 
lated with a ir phase shift; specially in the region just 
before v = 2. The results of Fig. 1 can help to interpret 
the experiments of Ref. Q and Ref . ||] , which observed 
anticorrelated behavior in the traces and spin alternation 
with increasing B, respectively. 

Having analyzed the non-interacting model we shall 
next estimate Coulomb interaction effects by including 
the selfconsistent Hartree potential 
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where k is the semiconductor dielectric constant (k = 
12.4 for GaAs). Electronic exchange and correlation en- 
ergies will be considered using the local-spin density ap- 
proximation (LSDA) within a spinorial formalism, since 
Hd breaks the symmetry of a single spin quantization 
axis. The LSDA relies on Monte Carlo calculations for 
the non polarized and fully polarized electron gases at 
B = |Q. A functional theory including current den- 
sity dependence, better adapted to systems in a mag- 
netic field, is known to exist [ fL5| . Nevertheless, current 
density corrections are quite small for moderate mag- 
netic fields and, besides, the numerical solution of the 
current-density-functional equations is almost unfeasible 
in a symmetry unrestricted case, unless strong smooth- 
ing approximations are introduced In the results 
shown below we consider weak vertical magnetic fields, 
as quantified by the condition v > 8, where the filling 
factor v is obtained from the dot central density pc as 
v = 2n£ 2 B pc (with £b — hc/eB the magnetic length). 

Defining the spin density matrix in terms of the Kohn- 
Sham spinors {ifi(v, 77), i = 1, . . . , N} as 
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the LSDA exchange-correlation functional E xc [p vv i] 
yields the following 2x2 potential matrix 
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Details on the evaluation of the functional derivatives in 
LSDA can be found in Ref. Jl|]. The resulting Kohn- 
Sham equations read 
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We have solved the set of Eqs. (|Tl]) by discretizing 
the xy plane in a uniform grid of points and applying 
an iterative scheme to reach full selfconsistency in Vr (r) 
and V xc , vv '(r). In some cases this procedure might get 
trapped in a local minimum. Therefore, several calcu- 
lations with different random initial conditions have to 
be used to ensure that the proper energy minimum is 
reached. The stability with the number of mesh points 
has also been checked. 

Figure 2 displays the density p(r) and spin magneti- 
zation m(r) for N — 9 electrons in a circular confining 
potential with ui x = ui y = 6 meV. We have selected this 
electron number as a representative case where to check 
the robustness of the analytically predicted spin inver- 
sions. A SO coupling with zq — 62 A and B = 2.5 T have 
been assumed. The magnetization density indicates the 
local orientation of the spin vector and it is related to the 
spin-density matrix by m x = 2Re[pfjJ, m y = 2Im[pjJ 
and m z — pff — p^. We note from Fig. 2 that circu- 
lar symmetry is conserved by both p and m Zl while the 
horizontal magnetization rri|| = (m x , m y ) shows an angu- 
lar dependent texture. This result is in good agreement 
with the analytical solution given above, that predicts 
m||(r) ~ p(r) (y,x). Note also that for this particular zq 
and B the vertical spin is predominantly inverted, giving 
a negative value for the total vertical spin (S z ). 

Figure 3 shows (S z ) as a function of B and the inten- 
sity of the SO coupling, given by zq, for the same dot 
of Fig. 2. In agreement with the above discussion, for 
B < 2.2 T the noninteracting model predicts spin in- 
version when decreasing the dot width. The LSDA also 
yields spin-inverted regions although with some conspicu- 
ous differences that can be attributed to a higher rigidity 
in the electronic structure. The interaction inhibits the 
spin flip at low magnetic fields and low widths shifting 
the inversion region to 1.8T < B < 2.6T and leaving 
only a small residue for zq < 40 A and B < 1.2 T. It is 
worth to mention that although in the laboratory frame 
(S z ) is not restricted to discrete values (because of the 
transformation C/i), in practice its fluctuations increase 
with the SO strength but they are generally small. In 
Fig. 3 the deviations from ±Ti/2 when zq = 62 and 48 
A are ~ 5 % and ss 20 %, respectively. Experimentally 
the dot width can be controlled with the electric gates; 
therefore, the results of Fig. 2 suggest a spin switch be- 
havior, controllable with vertical magnetic and electric 
fields. 

In Fig. 4 we show the vertical spin as a function of 
the applied magnetic field and the dot deformation, for 
two different values of the SO coupling. In this figure the 
mean value (uj x + w y )/2 was kept fixed to 6 meV while 
the ratio 5 = lo v /u> x was varied to obtain different el- 
liptical shapes. Comparing left and right panels we see 
again an interaction-induced quenching of the spin in- 
version at low magnetic fields. Although this occurs for 
the two displayed widths, at zq = 48 A a larger region 



with inverted spin is found. Figure 4 shows that, hav- 
ing fixed the SO coupling strength, spin inversion can be 
achieved in many cases either by increasing the deforma- 
tion or, alternatively, by increasing the magnetic field. 
The comparison of the N — 9 numerical results with the 
analytical model allows us to conclude that, in spite of 
the differences, the spin inversions in LSDA are quali- 
tatively similar to the analytical ones. As a final piece 
of information we mention that the energy gap between 
the highest occupied and lowest unoccupied Kohn-Sham 
leves in the above cases stays in the range [0.3, 0.9] meV. 
This result provides a mesure of the ground state sta- 
bility and, therefore, of the relative spin stiffness against 
thermal fluctuations. 

In summary, the mechanism of spin switch through 
SO coupling has been analyzed. It has been shown that 
SO coupling can lead to anticorrelated behavior in the 
B evolution of neighboring levels and to up and down 
spin oscillations; qualitatively similar to experimental ob- 
servations. The combined effects of SO coupling, weak 
magnetic fields and deformation have been studied; first 
analytically with a non interacting model and; second, by 
taking into account interaction effects within the LSDA. 
The diagrams with the spin dependence on these param- 
eters suggest the characterization of the quantum dots as 
spin switches, of relevance to spintronic technology. 
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FIG. 1. Single particle energies in modified atomic units 
(for GaAs 1H* « 12 meV) in the noninteracting model. Cir- 
cular symmetry with u) x = tu y — 1.1 meV has been assumed. 
Each line has been shifted vertically by a small amount rep- 
resenting the charging energy. Light and dark gray color cor- 
respond to up and down spin, respectively. 

FIG. 2. Density p and magnetization m for a circular dot 
with N = 9 electrons, Tiu> x — hu> y — 6 meV, B — 2.5 T, and 
«o = 62 A. The values have been scaled by 1/9 of the dot 
central density, i.e., p = (2.2 10~ 5 A~ 2 )/9. 

FIG. 3. Vertical spin evolution in the B-zo plane for the 
same quantum dot of Fig. 2. White (black) color indicates 
upward (downward) total spin. 

FIG. 4. Vertical spin evolution in the B-S plane, where 
S = LUy/cux labels the deformation, for a fixed value of the SO 
coupling. Upper row corresponds to zq — 62 A in the non 
interacting model (left) and the LSDA (right). Lower row 
shows the same results for z = 48 A. We have used the same 
color convention of Fig. 3. 
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